Dynamical model for the formation of patterned deposits at receding contact lines 
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We describe the formation of deposition patterns that are observed in many different experiments where 
a three-phase contact line of a volatile nanoparticle suspension or polymer solution recedes. A dynamical 
model based on a long-wave approximation predicts the deposition of irregular and regular line patterns 
due to self-organised pinning-depinning cycles corresponding to a stick-slip motion of the contact line. 
We analyze how the line pattern properties depend on the evaporation rate and solute concentration. 

PACS numbers: 68.15.+e, 47.57.-s, 81.15.Lm, 81.16.Rf 



The last decade has seen huge growth in interest in phe- 
nomena that accompany evaporative and convective dewet- 
ting of suspensions and solutions. Well known are the de- 
tailed studies of the coffee stain effect [Ql 01 that analyse 
the deposition and resulting structures left behind by a re- 
ceding three-phase contact line of an evaporating drop of 
suspension upon a solid substrate. In particular, Ref. J2l 
describes a large range of different deposition patterns in- 
cluding cellular and lamellar structures, single and multi- 
ple rings, and Sierpinski gaskets. Other observed struc- 
tures include crack J3|] and chevron |0] patterns. Recently 
it has been shown that evaporating polymer solutions ]5|- 
7] and (nano)particle suspensions [f8l-ll0h may be used to 
fabricate strikingly regular stripe patterns, where the de- 
posited stripes are parallel to the receding contact line 
and have typical distances ranging from 10-100/im. The 
goal is to use this effect as a non-lithographic technique 
for covering large areas with regular arrays of small-scale 
structures, such as, e.g., concentric gold rings with poten- 
tial uses as resonators in advanced optical communications 
systems 11 ill . The deposited patterns from more complex 
fluids, such as polymer mixtures Il2ll and DNA solutions 
1 13ll . are also investigated. The occurrence of regular stripe 
patterns is a somewhat generic phenomenon, that is not 
only observed for different combinations of substances but 
also in a variety of experimental setups that allow for slow 
evaporation. Examples include the meniscus technique in 
a sphere-on-fiat geometry 0, 01. a controlled continuous 
supply of liquid between two sliding plates to maintain a 
meniscus-like surface [Q] and dewetting forced by a pres- 
sure gradient ifioll . Interestingly, besides the stripes paral- 
lel to the receding contact line, a variety of other patterns 
are observed, including regular orthogonal stripes |@], su- 
perpositions of orthogonal and parallel stripes yfl, regular 
arrays of drops JH [Hfl and irregularly branched structures 
lfl4lTl5ll . This behaviour is highly sensitive to the particular 
experimental setup and parameters. 

Despite the extensive number and variety of experi- 
ments, an explanation of the formation of the regular pat- 
terns has been rather elusive. Although most studies agree 
that the patterns result from a stick-slip motion of the con- 
tact line caused by pinning/depinning events iBr&lfll, 161 
no dynamical model of the periodic deposition process ex- 
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FIG. 1 : Sketch of a liquid front that recedes due to evaporation 
and convection with a varying velocity v(t). The deposition pro- 
cess is characterized by the film thickness profile h(x, t), the con- 
centration profile 4>(x,t) (in the bulk film), and the particle layer 
thickness h p (x, t) — <j>{x, t)h(x, t) (outside the bulk film). 



ists. Most models assume a permanently pinned contact 
line (see 017l Il8ll and references therein) and are therefore 
only able to describe the formation of a single line deposit. 
A non-isothermal Navier-Stokes simulation shows depin- 
ning from such a single line but no periodic deposits lllql . 
The model of Ref. H% describes drop arrays formed via di- 
rected dewetting of the solvent that are subsequently dried. 

In this Letter we discuss a generic close-to-equilibrium 
model for the evaporative and convective receding of a 
three-phase contact line of a solution or suspension on a 
solid substrate. We show that solely having a viscosity 
that diverges at a critical solute concentration is sufficient 
to trigger a self-organised periodic pinning-depinning pro- 
cess that results in the deposition of regular line patterns. 
The model can easily be extended to incorporate other pro- 
cesses and in the future may be employed to assess their 
influence on the basic mechanism that we describe here. 

We consider a thin film of an evaporating partially wet- 
ting nanoparticle suspension (or polymer solution) in con- 
tact with its vapour on a flat solid substrate (see Fig.[Ql. As- 
suming that all surface slo pes are small, one may employ a 
long- wave approximation EOll and derive two coupled evo- 
lution equations for the film thickness profile h(x,t) and 
the vertically averaged solute concentration field (f)(x, t): 



ft, 



d t h = d x [Q(h,<f>)d x p(h)] - -ip(h) - tip), (1) 

P 

d t {ct>h) = d x [4>Q{h, 4>)d xP {h)] + 3 X [D{(f>)hd x (t>] (2) 



The mobility Q(h, <j>) = h 3 /3r]((j)) models Poiseuille flow 
and incorporates the dynamic viscosity that exhibits 
a strong non-linear dependence on the local solute concen- 
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tration. We employ the Krieger- Dougherty law 1211122 

where r/ is the viscosity of the pure solvent. We have 
scaled (f) by the concentration at random close packing 
((j> c = 0.63) so that r\ diverges when 0—^1. The pre- 
cise value of the exponent v = [rj\ <p c depends on the type 
of suspension. For non-interacting particles (i.e. particles 
that have no net attractive forces between them and only 
have excluded volume interactions), values for v between 
1.4 and 3 are discussed, depending on the shape of the par- 
ticles 12 in . For spherical particles the factor [77] = 2.5, 
givi ng y = 1.575. Other thin film models use v = 2 
1 19l 231. For interacting solute particles, values for v as 
low as 0.13 are reported [24]. Depending on the particular 
system, the transition at (j) c is either referred to as jamming 
or gelation | f24} |. Here we fix v = 1.575, although we have 
found that the effects we describe below are even stronger 
for smaller v. 

The first term on the right hand side of Eq. (TJ (con- 
served part) corresponds to convective transport of the liq- 
uid whereas the second term (non-conserved part) models 
evaporation. The convective flow is driven by the gradient 
of the pressure 



p(h) = --yd xx h - 11(h), 



(4) 



where the first term is the Laplace pressure (7 is the surface 
tension) and the second is the disjoining pressure 11(h) = 
2S LW d 2 /h 3 +S p exri-(h-d )/l ]/l that models a par- 
tially wetting fluid l!25L l26ll . Here, l is the Debye length, 
d is a molecular interaction length, S LW = —AjYlixd 2 ^ 
and S p < are the apolar and polar spreading coeffi- 
cient, respectively, and A < is the Hamaker constant. 
To derive the second term in Eq. ([TJ we assume the sys- 
tem is close to equilibrium and near to saturation and so 
evaporation is slow. In this limit evaporation with a rate /3 
is driven by the difference of the scaled pressure pip and 
the chemical potential of the ambient vapour p 11271 |28h . 
Latent heat effects may be neglected, and the density p 
is assumed to be equal for particles and solvent. The 
first and second terms on the right hand side of Eq. © 
model convective and diffusive transport of the particles, 
respectively. Note that the diffusion coefficient depends on 
concentration and we employ the Einstein-Stokes relation 
D(4>) = kBT/6irr 7]((j)), where ks is the Boltzmann con- 
stant, T the temperature, and r the particle radius. 

Models related to Eqs.JT) and (O are used in studies of 
particle-laden film flow [23fl (without evaporation or wet- 
tability effects) and dewetting of suspensions of surface 
active particles fl9i] (different wettability regime). In the 
limit <p — > 0, our theory reduces to that used in ]27|] to 
study the fingering instability of an evaporative front of a 
pure liquid. We choose our scaling and some of the param- 
eters to be the same as in rf27ll : The dimensionless chemi- 
cal potential M = pp/\S p \ = —0.003 and the diffusion 
number D = 3k B T/r [67rA 2 \S p \] 1 / 3 = 0.0003, where 




5000 10000 15000 





| 1 1 1 1 1 1 u 
1 fl>) ! 








70000 75000 80000 



35000 40000 45000 50000 55000 
X 



FIG. 2: (color online) Typical deposit profiles: (a) a single line, 
(b) a finite sequence of lines, (c) an intermittent line pattern, and 
(d) a regular line pattern. The parameters are <f>o = 0.41 and from 
(a)to(d): Qq = (14.7,0.1,0.147,0.464) x 10~ 6 . 



S p = S p exp(d /l )/lo- Our main control parameters 
are the initial mean (dimensionless) concentration (j> and 
the evaporation number f2 = 187r/3?7o7//o[67ry4 2 |S' p |] 1 / 3 
which represents the ratio of the time scales for convec- 
tion and evaporation of a film without solute. We solve the 
nondimensional model by discretizing over a spatial do- 
main of finite length L. Deposited patterns are obtained by 
direct time simulations using a variable-step variable-order 
backward difference scheme starting from an initially step- 
like front that becomes smooth in the early evolution 12911 . 

In the evaporative dewetting process one encounters 
different types of receding fronts. In the case with- 
out solute 12711 one may distinguish the limiting cases of 
(i) convection-dominated and (ii) evaporation-dominated 
dewetting for small and large values of Q , respectively. 
In case (i) the front recedes rapidly and convective motion 
maintains a capillary ridge despite evaporation. In case (ii) 
convection is much slower than evaporation, the front re- 
cedes slowly and there is no capillary ridge. In the pres- 
ence of a solute the situation is more complex and station- 
ary receding front shapes are not found for some parameter 
values. In general, on starting with a uniform concentration 
4>(x, 0) = (fto, as the front recedes, it deposits part of the 
solute in a smooth layer. Evaporation in the contact line re- 
gion increases the local concentration (f> and consequently 
also the viscosity. When (f) — )■ 1, the convective motion 
in the contact region stops completely (the suspension be- 
comes locally jammed) due to the strong non-linearity in 
Eq. (O. One may say that the initial convection-dominated 
stage changes into an evaporation-dominated regime as the 
local evaporation number Q = (rj((fi)/r] )no becomes 
large. At this stage the front seems pinned. However, it 
is actually still moving extremely slowly, solely by evapo- 
ration, and deposits a first line of solute. During this, the 
local concentration in the contact region decreases, i.e., Q 
decreases, until the front depins and convective motion re- 
sumes. Subsequently, after the initial line is deposited, var- 
ious scenarios are possible: (a) a deposit of uniform thick- 
ness (i.e. only a single line is deposited), (b) deposition 
of a finite number of lines followed by a layer of constant 
thickness, (c) intermittent line pattern, and (d) regular line 
pattern. In the narrow region of parameter space where the 
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FIG. 3: (color online) Morphological phase diagram of deposi- 
tion patterns, in the plane spanned by evaporation number fl$ 
and bulk concentration 4>o- Symbols denote performed simula- 
tions. In the central region we find (d) regular line patterns (o) 
and outside of this region (*) we observe: (a) single lines, (b) 
multiple lines, (c) intermittent patterns, and (e) no lines. For typ- 
ical deposit profiles for (a)-(d) see Fig. [2] Results from along the 
dashed lines are analysed at Figs.|4]and|5] 



'limiting' case (c) occurs, we find that this behaviour is nor- 
mally very sensitive to computational details (which is not 
the case for (d)), so we believe that the intermittent line pat- 
terns represent a 'chaotic deposition'. Typical profiles are 
displayed in Fig.|2]and the corresponding parameter ranges 
are marked in the phase diagram Fig. [3j 

Ignoring initial transients, we distinguish two main types 
of deposits: those of uniform thickness (regions (a), (b) and 
(e) in Fig. [3) an d periodic line patterns (region (d)) . In the 
latter one observes (as in the experiments |@, [[(J, [111]) a 
regular stick-slip motion of the contact line, since the typi- 
cal time scales for convection-dominated and evaporation- 
dominated front motion may differ by orders of magnitude. 

Structure formation results from a subtle interplay be- 
tween convection, evaporation, and diffusion. The basic 
mechanism of line deposition described above stems from 
a balance of convective and evaporative motion. Diffusion 
does not change that picture as long as its time scale is 
not much shorter than the convective and evaporative time 
scales. Increasing D Q merely shrinks the region (d) in 
Fig. [3] However, in the unlikely situation that diffusion is 
so fast that nanoparticles diffuse away ahead of the reced- 
ing contact line (e.g. when D = 0.3), then line deposition 
is suppressed. We do not consider this case here. 

Next, we analyse re gula r line patterns as obtained from 
long-time simulations 112911 - Excluding the initial transient 
we take a sequence of N regular deposition periods (lines), 
where 10 < N < 100, depending on the period of the de- 
posit and required CPU time. We measure the amplitude, 
relative width (defined as 2a / [period] , where a is the stan- 
dard deviation) and skewness of the lines, and the period 
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FIG. 4: (color online) Left: dependence of the regular line pattern 
properties on (f>o in the range bounded by the vertical dashed lines 
that corresponds to region (d) of Fig. [3] Q$ = 4.64xl0~ 7 is 
fixed. Right: Line patterns for <j> = 0.305, 0.32, 0.41, 0.475 and 
0.499 (from the top) indicated by dotted lines in the upper left 
panel. 



of the line pattern. We find that these quantities strongly 
depend on both the evaporation number f2 an d the con- 
centration (p . We focus on two cuts through region (d) in 
Fig. [3] (dashed lines). 

First, we fix f2 = 4.64 x 10 -7 and vary the bulk con- 
centration (j) . Fig. |4] presents line characteristics and se- 
lected profiles. On increasing (fr from a region without pe- 
riodic line deposition, one first finds large amplitude almost 
solitary peaks separated by very large distances. On fur- 
ther increasing (j) , the amplitude first hardly changes and 
later decreases. The period continuously decreases while 
the relative width increases. For higher <J) the deposit pat- 
tern becomes almost uniform, with a small amplitude har- 
monic modulation. Finally, the amplitude goes to zero (at 
finite period) at the upper border of the region (d). 

Second, Fig. [5] presents results for fixed <^> = 0.41 and 
varying i? . Increasing f} , moving from the narrow re- 
gion (b) of multiple lines (Fig. [3]), one passes through a 
very narrow band of intermittent line patterns (c) followed 
by the region (d) of regular line patterns. For the lowest 
values of f2 in (d), the patterns have a relatively small 
period and a small but non-zero amplitude. The strongly 
anharmonic peaks are skewed to the right with their tail 
pointing towards the receding front. On increasing f2 the 
period increases. The amplitude, however, first increases 
and then decreases, until at a certain threshold the pattern 
ceases to be periodic and we arrive in the narrow border re- 
gion (b), where only a finite number of lines are deposited. 
Correspondingly, the relative width decreases, as lines get 
more peaked. Remarkably, the skewness changes sign, i.e. 
the tail of the lines shifts from pointing towards the reced- 
ing front, to pointing away. This effect was observed in 
experiments on nanoparticle suspensions l3fj|| and can be 
explained as follows: For smaller values of f} , the cap- 
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FIG. 5: (color online) Left: dependence of the regular line pattern 
properties on J? in the range bounded by the vertical dashed lines 
(region (d) of Fig. [3}. 4>o = 0.41 is fixed. Right: Line patterns for 
J?o = (0.17, 1.0, 2.15, 5.27, 5.66) x 10~ 6 (from the top) indicated 
by dotted lines in the upper left panel. 



illary ridge is large and accommodates a large amount of 
nanoparticles as it recedes. When the front pins, the cap- 
illary ridge is evaporated and these nanoparticles are de- 
posited in the thick tail to the right. When the liquid front 
depins, it results in a further drop in the deposition thick- 
ness (seen as shoulder in the tail). For higher values of f2 
the capillary ridge is smaller and so the right tail is smaller. 

In conclusion, we have studied a generic model for 
the close-to-equilibrium deposition dynamics onto a sur- 
face from a polymer solution or nanoparticle suspension. 
The model incorporates wettability, capillarity, evapora- 
tion, convective transport of the solution and diffusion of 
the solute and has been derived employing a long-wave ap- 
proximation. We find that a strong nonlinear dependence 
of viscosity (i.e., the front mobility) on concentration trig- 
gers, in an intricate interaction with evaporation and dif- 
fusion, the deposition of periodic and aperiodic line pat- 
terns as observed in experiments for many different mate- 
rials and settings lBl5l 4TTl[30h . We believe that the model 
explains a basic mechanism for the formation of regular 
line patterns. They result from a self-organised cycle of 
deposition-caused pinning-depinning events that is experi- 
mentally often described as a stick-slip motion fl flollTlll . 
In the future the basic dynamical model should be extended 
to study the influence of other important effects on the de- 
position. These include thermal effects, solute-dependent 
wettability and the nature of the solvent-solute interaction. 
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